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We derive the stochastic master equations, that is to say, quantum filters, and master equations 
for an arbitrary quantum system probed by a continuous-mode bosonic input field in two types 
of non-classical states. Specifically, we consider the cases where the state of the input field is a 
superposition or combination of: (1) a continuous-mode single photon wave packet and vacuum, 
and (2) any number of continuous-mode coherent states. 
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I. BACKGROUND AND MOTIVATION 

The production and verification of non-classical states 
of light, such as single-photon states [T] and superposi- 
tions of coherent states (also known as Schrodinger cat 
states) [IHlj, has become routine. In particular, the pro- 
duction of single photon states has been achieved in a va- 
riety of experimental architectures such as: cavity quan- 
tum electrodynamics (QED) [51[3], quantum dots in semi- 
conductors [7], and recently in circuit QED [8 . Such 
non-classical states have been considered in connection 
with quantum computing [21 llOj and secure communica- 
tion [TT] over quantum networks [T^]. 

A basic problem in quantum optics concerns the ex- 
traction of information about a system of interest (two- 
level atom, cavity mode, etc) from light scattered by the 
system. Figure [11 Based on measurements of the scat- 
tered, or output, light, one can determine a conditional 
state from which one can make estimates of observables 
of the system. A general approach to estimation prob- 
lems of this kind, called /iZtermi? problems, was developed 
by Belavkin J13j-[16j within a framework of continuous 
non-demolition quantum measurement in the case where 
the input probe field, B{t) in Figure R] is a quantum 
white noise with vacuum state (or more generally Gaus- 
sian state, see [ITj- [2Q]). Belavkin's formulation, which 
generalizes the classical nonlinear filtering theory [3T], 
is quite general. For example, in the schematic repre- 
sentation of a continuous measurement process shown in 
Figurelll the measurement signal Y{t) produced by a de- 
tector (e.g. photon counter or homodyne detector) may 
be the number of quanta in the output field, or alterna- 
tively it may be a quadrature of the output field. 

One obtains a filtering equation which is a stochas- 
tic differential equation of the state Q(t) conditioned on 
Y{t): in the later terminology employed in quantum op- 
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FIG. 1. A schematic representation of a continuous measure- 
ment process, where the measurement signal produced by a 
detector is filtered to produce estimates X{t) = nt{X) = 
tr[p{t)X] of system operators X at time t. 



tics, the output is referred to as a quantum trajectory 
and the filtering equation as a stochastic master equation 
P^24| . Averaging over the measured output is equiva- 
lent to a non-selective measurement, and the correspond- 
ing state will satisfy the corresponding master equation. 
The choice of detection scheme on the output field deter- 
mines the particular selective evolution, usually referred 
to as an unravelling of the master equation in quantum 
optics. To date, quantum trajectories and quantum fil- 
tering have only been developed for input fields that are 
in a Gaussian state, with specific cases being coherent 
state fields (this includes vacuum fields as special case), 
thermal fields, and squeezed fields [25142 7j . While the re- 
sulting equations allow us to estimate non-commutating 
observables of the monitored system, the Gaussian nature 
of the inputs ensure that they appear formally similar to 
the classical equations. The aim of the present paper is 
to extend the theory to classes of non-classical inputs. 

In this article we extend Belavkin's quantum filtering 
theory and the input-output theory of quantum optics 
[22l [27] to non-Gaussian continuous-mode states pudd 
which are superpositions or combinations of 

a) a continuous-mode single photon and vacuum, and 

b) continuous-mode coherent states, i.e. continuous- 
mode cat-states. 

The problem to be tackled here is to derive the master 
and stochastic master equations for a "system" (g) "field" 
with initial state po ^Pficid and unitary evolution process 



U (t) when the field state is one of the states above. To 
make the problem tractable, we seek a larger representa- 
tion of the form "extended system" (g) "field" where 

"extended system" = "ancilla" (E) "system" 

such that, for all system observables X, 

trsystem(g.fiold JPO «> Pfiold U (t)^ {X ® I) U {t)\ 
^ trancilla(g)system(g)field\Pa ® Po ® -i vac 

U(t)\R{t)®X®I)U{t)} (1) 

where U (t) is a unitary evolution process coupling the 
ancilla, system, and field, pa is a fixed state of the an- 
cilla, Pvac — |0)(0| is the vacuum state (projection) for 
the continuous-mode field, and R (t) is some process tak- 
ing values in the observables of the ancilla. The filter- 
ing problem may then be solved for the extended system 
with reference to the vacuum state for the field using 
traditional techniques. 

The extension to single photon states is interesting for 
foundational reasons [25] as well as the aforementioned 
technological reasons [S]. Likewise, quantum filtering for 
cat states is of foundational importance, while practical 
uses would be towards quantum enhanced metrology [29] . 
One possible application would be to quantum enhanced 
metrology of a time varying parameter [301 131] . 

This article is structured as follows. In Section |ll] 
we review standard input-output theory. Specifically 
we consider the idealized quantum white-noise model of 
a quantum stochastic differential equation (QSDE) and 
use it to derive the master equations and quantum tra- 
jectories for Gaussian fields. Then we review a general 
parametrization to specify the system environment cou- 
pling for input-output systems. Using this parametriza- 
tion we review the methods, recently introduced [32ti34j . 
to simplify and formalize the network theory of cascaded 
open quantum systems and quantum feedback networks. 



Section III is focused on deriving the master equation 
and stochastic master equation (quantum filter) driven 
by continuous-mode single photon wave packets. We 
generalize the single photon filter to any superposition 
or combination of single photon and vacuum input field. 
The system that is probed is left arbitrary so in general 
our filter can apply to qubits, qudits and mechanical os- 
cillators. As an example we calculate the single photon 
filter for a two level atom (or qubit) dispersively cou- 
pled to the field. We derive the trajectories for both 
a homodyne type measurement and a photon counting 
measurement. 



In Section IV we present the extension to superposi- 
tions of coherent states. We derive the cat-state-filter 
for an arbitrary quantum system and an arbitrary cat 
state. Again we illustrate the filtering equations with a 
qubit system and homodyne and photon counting mea- 
surements. 

In Section [V] we conclude and discuss our future re- 
search and some open questions. 



Notation The commutator and anti-commutator will be 
denoted as [A, B]= AB- BA and [A, B]+ = AB + BA, 
respectively. We set VaB = A^BA~ i(^t^S + BA'fA) 
and V\B = ABA't - ^{A^ AB + BA^^A). 

The scattering, coupling and Hamiltonian operators 
describing a given Markovian open system coupling will 
be written as a triple G = {S,L,H), to be explained in 
more detail in Section [lI A[ and this provides an operator- 
valued parameterization of the system. The associated 
superoperators are 

Lindbladian : CgX = —i\X, B\ + VlX, 
Liouvillian : CqP = —i[H, p] + "D^p, 

and note that, for traceclass p and bounded X, 

tiipCaX} = tv{XC*aP}. 



II. MODELS OF OPEN QUANTUM SYSTEMS 

In this section we briefiy review quantum stochas- 
tic calculus (input-output theory) and quantum filtering 
(trajectories) for a system coupled to a heat bath mod- 
elled as a boson field in the vacuum state. 



A. Input-Output Model Using QSDEs 



Hudson and Parthasarathy [351 136] showed how to di- 
late a dissipated completely positive semigroup evolution, 
with Lindblad generator, to a unitary model on the sys- 
tem space with a (Bose) Fock space ancilla. Here they 
developed an analogue to the Ito theory of stochastic in- 
tegration with respect to creation, annihilation and scat- 
tering process B^ (t) , B (t) and A (i) . They showed the 
existence and uniqueness of solutions to unitary quantum 
stochastic differential equations of the form 



dU{t) =<{ (5 - 1) dA (t) + LdB\t) 

L^SdB{t) - {\L'^L + iH)dt\u{t). (2) 



where 



G={S,L,H) 



consists of a unitary S describing photon scattering 
phase, a bounded operator L describing coupling to the 
creation mode of the field, and a bounded Hermitean 
operator H describing the system Hamiltonian. (The 
result has been extended to non-bounded coefficients.) 
The increments are future pointing operator-valued Ito 
increment, that is dB{t) = B{t + dt) — B{t) is a for- 
ward of the quantum noise. In particular, we have. 



[U{t),dB{t)] = [U{t),dB^{t)] = [U{t),dA{t)] = 0. The 
full quantum Ito table is 



X \ dt dB dA dB'f 
dt Q Q 
dB dB dt 
dA dA dB'' 
dB^ 



(3) 



More generally, for quantum stochastic integral processes 
X(t), Y{t), one has the Ito product rule 

diX{t)Y{t)) = idX{t)) Y{t) + X{t) dY{t) + dX{t) dY{t). 

Independently, Gardiner and CoUett developed an 
equivalent quantum input-output theory [26l [27] based 
on Lehmann-Symanzik-Zimmermann scattering theory 
of Bose white noise processes. Formally one begins with 
singular fields satisfying 

with the connection to the regular processes being for- 
mally 

si' [t) = / &^ (s) ds, B{t)= f b (s) ds, 



A(i) = / b^s)b{s)ds. 
Jo 

The quantum stochastic calculus may then be under- 
stood as effectively arising through Wick ordering of the 
singular fields. 

The multiple input version is relatively straightfor- 
ward. We have n independent inputs bj and with Bj (t) = 

/o ^3 (^) ^^' ^Jk (^) — /o K (■*) ^fc (^) '^^^ 6^*^-' ^^ have 



where we now have parameterizing operators 



S = 



/Si 



\ Snl 



Sir, 



On 



,L = 



fLi 



\L„ 



,H 



with 5' unitary and H self-adjoint. For simplicity we 
treat the case of a single input and output. 



B. Heisenberg-Langevin Equations 

The Heisenberg dynamics of arbitrary system operator 
X is defined by transforming to the Heisenberg picture 

jt{x)^uHt){x®hM)u{t). 



(We will usually drop the subscripts "system" and "field' 
when there is no confusion.) From the quantum Ito prod- 
uct rule and table one deduces the QSDE for a system 
operator jt{X) — X(t): with all system operators trans- 
formed to the Heisenberg picture. 

djtiX) ^ jt{CX)dt 

+ dB\t)jt{S^X,L])+jt{[L^,X]S)dB{t) 
+ jt{S^XS-X)dAit). (4) 

C. Derivation of the Master Equation 

Suppose that the system is in an initial state p(0) — po 
and that the joint state of the system and bath is /9o®-Pvac 
where Pvac = |0)(0| is projection onto the vacuum state 
of the field. The state of the system, g{t), obtained by 
averaging over the environment at a given time t is then 



g{t) - trfield [Uit){po ® Pvac)C/^(i)] 



(5) 



We wish to obtain a differential equation for the average 
of an observable X of the system at time t: 

^i {X) = trgystcmigjfiold {jt {X) Qo <8) Pvac} 
= trsystcm {P (t) X} , 

and from the Heisenberg-Langevin equation Q we have 

dwt (X) = trsystcmigifiold {djt {X) Po (K) Pvac} 

— trsystem«)field {jt {^gX) Po (Xl Pvac} dt, 

as the increments dB,dB',dA vanish in the vacuum 
state. We therefore obtain the equation 

dvjt (X) 



dt 



= ^t i^cX) , t^o {X) = trgystcm {pqX} 



which may then be expressed as the master equation 
dg(t) 



dt 



Chgit) = -i[H,git)]+Vlg{t), (6) 



with initial data po- Note that the master equation m 
is a consequence of the QSDE model. 



D. The Input-Output Relations 

The output field -Bout is obtained from the input by 
moving into the Heisenberg picture: 

Bout (t) = U (i)^ (/.ystcm ® B (i)) U {t) 
= C/(r)^(/,y,tem®5(i))f/(T) 

for any t > t. Again from the quantum Ito calculus we 
find 



dBout (t) ^ Jt (S) dB (t) + Jt (L) dt. 



(7) 



Note that the output field again satisfies the canonical 
commutation relations. 
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FIG. 2. An open quantum system. The input field (before 
interaction) is represented by the operator B(t) and output 
field (after interaction) is denoted by -Bout(i). 



E. Derivation of the Quantum Filter (Stochastic 
Master Equation) - Quadrature Case 

We suppose that we continuously monitor the quadra- 
ture phase using perfect (100% efhciency) honiodyne de- 
tection. This entails measurement, for each t > 0, of the 
field 

r(t)-i?out(t)+SLW 

= [/(i)^ (/system ®Q(i))C/(t) 

where Q{t) = B (t) + B^ {t). We note that the set of ob- 
servables {Y (t) : t > 0} is self-commuting and we may si- 
multaneously diagonalize (and measure!) all observables. 
At any time t, we may additionally estimate an observ- 
able that commutes with the observables up to time t. 
This includes observables X (r) for t > t, since 

[XiT),Y{t)] 

= U{T)*[X(g) /field, /system Q {t)] U (t) 

eeO. 

This is the non-demolition property. Quantum filtering 
is the estimation of jt {X) based on observations of the 
output processesjy (s) : < s < i}. Fig. fl] depicts the 



scenario we are considering. From the Ito calculus we 
see that 

dY{t) = {L{t) + L\t))dt + dQ{t). 

Defining the expectation 

£[•] =tr{po«)Pficid(-)} 
for a given state po ® Pficid, we seek to minimize 

n{x{t)-3t{xm 

over all observables X (i) in the algebra J^j generated by 
{Y (s) : < s < i}. The minimizer is called the least- 
squares estimator for X (t) given {Y (s) : < s < t} and 
will be denoted as 



X{t)^7rt{X)=E[UX)\yt] 



(8) 



The later notation suggest that in tt^ (X) is the condi- 
tional expectation oi jt {X) given the past history, which 
would be the classical interpretation. While conditional 
expectations generally do not exist in the quantum prob- 
abilistic setting, the nondemolition property above suf- 
fices to allow one to realize precisely this interpretation, 
see for instance [HI |37] . The conditional expectation can 
indeed be interpreted as an orthogonal projection onto a 
subspace of commuting operators yt- This means that 
jtiX) — 7rt(X) is orthogonal to this measurement sub- 
space yt, that is. 



E[ijt{X) - 7rt{X))C] ^ 



(9) 



for all operators C belonging to the measurement sub- 
space yt, [TU]. Setting C ^ I shows that 

E[7Tt{X))]=E[jtiX)]. 



Quantum Stochastic Differential Equation 

dU{t) = {LdB^ - L^SdB - {^L^L + iH)dt]U{t). 

dX{t) =/:t[X]dt + dB^[X,L] + [L\X]dB 

dV = {L + V)dt + dB + dB^ 



Quantum expectation 



Quantum conditional expectation 



IVlaster Equation 

dg= -idt[H,g\ + dtV^[L]g 



classical 
expectation 



dp -■ 



Stochastic Master Equation (Quantum Filter) 

-I dt[H, p] + dt Pt [L\p + {dV - tr[(L + L'^)p])dtn[L]p 



FIG. 3. Relationship between QSDEs for the unitary, system and environment; master equation; and stochastic master 
equation (filter) for open quantum systems ("Hl is given by (14 1. The difference between the master equation and stochastic 



master equation is due to the difference in the type of expectation take, i.e. unconditioned or conditioned, respectively. 



Now let us return to the vacuum state for the field: 



Pfiold 



We shall recall a simple derivation of the 



filter using an analogue of a the characteristic function 
technique of classical filtering jSS] . We introducing a pro- 
cess C (t) satisfying the QSDE 



dC (t) = g (t) C (t) dY (t) , 



(10) 



with initial condition C (0) = /. Here we assume that g 
is integrable, but otherwise arbitrary. The technique is 
to make an ansatz of the form 



dnt (X) = atdt + Pt {X) dY (t) 



(11) 



where we assume that the processes at and /3t are 
adapted and lie in J^^. These coefficients may be deduced 
from the identity 

E[{nt{X)~jt{X))C{t)]^0 

which is valid since C{t) is in yt- We note that the Ito 
product rule implies I + II + III = where 



I = E[{d7Tt{X)-djt{X))C{t)], 
= E[atCit)+l3tjt{L + L^)Cit)]dt 
-E[jt{CGX)Cit)]dt, 
II = E[{7rtiX)-jtiX))dC{t)], 

= E [{TTt {X) - jtiXj) g (t) C {t)jt {L + it)] dt, 
III = E [{diTt {X) - djt {Xj) dC (t)] 

= E [Ptg [t) C {t)] dt + E [g{t)jt{[L\X])C{t)\ dt. 

Now from the identity I + II + III = we may extract 
separately the coefficients of g (t) C (i) and C (i) as g (t) 
was arbitrary to deduce 

TTt {{TTt {X) - 3t{X)) Jt (L + Lt)) + TTt Wt) = 0, 

TTt {at + PtJt {L + L^) - Jt (CgX)) = 0. 

Using the projective property of the conditional expecta- 
tion T:t{{TTtX)) = TTtiX) and the assumption that at and 
l3t lie in yt, we find after a little algebra that 

pt = nt {XL + L^X) - ^t {X) ^t {L + L^) , 
at = nt [CgX) - Pt^^t {L + L^) , 



so that the equation ( 11 ) reads as 
dnt{X) ^ nt{.CGX)dt 



(12) 



+{Tit{XL + L^X) -Trt{L + L^)TTt{X))dW{t), 

where the innovations process W(t) is a Wiener process. 
It is related to the measurement process Y{t) by the equa- 
tion 



and has the interpretation as given the difference between 
the observed change dY{t) and the expected change 
TTt{L + L^)dt in the measured field immediately after time 
t. Note that the increment dW{t) of the innovations pro- 
cess is independent of Trs{X) for all < s < i. 

It important to note that Q{t) (equivalent to a Wiener 
process) and W{t) (also a Wiener process) are distinct, 
and that Q(t) is not in the commutative observation sub- 
space. Some care is needed in interpreting equation ( 12 ) 
for the quantum filter. All of the terms in this equa- 
tion belong to the commutative subspace yt, and so (by 
the spectral theorem [TH]) are statistically equivalent to 
classical stochastic processes. 

The stochastic master equation may be expressed in 
terms of the density operator-valued stochastic process 

Pity. 

dp{t) = C*apit)dt + HLpit) dW{t). 
where we introduce 



UlP = Lp + L^p - tr { (L + L^) p] p. 



(14) 



The increments dW{t) can be generated independently of 
p{s) for all < s < t and the stochastic master equation 
above driven by the generated increments can thus act as 
a simulated quantum, trajectory of the state conditioned 
upon the measurement outcomes {j/(s); < s < t}. 



F. Photon Counting Case 

If instead we measure the number observable Y{t) = 

W{t)A{t)U{t) = Ao„t(i) = /g &L(s)6o,t(s)ds then the 
quantum filter is (see the survey paper [19| for the deriva- 
tion) , 



dp{t) = Cap[t)dt + JLp{t) dN{t) 



where 



JlP^ 



LpL^ 
t^ipUL} 



dY{t) = 7rt(i + L^)dt + dW{t) 



(13) 



and the innovations process in this case is given by 
dN{t) = dY — ti { p{t) L^ L}dt and is a compensated Pois- 
son process of intensity tT:{p{t)L^ L} . 



G. Cascade Connections 

A simple quantum network may be formed by con- 
necting the output of one system to the input of another 
system, (32l |34l |39l S^ • Fig. [I] illustrates the open quan- 
tum system G — {S,L,II) equivalent to the cascade of 
systems Gi = {Si,Li,Hi) and G2 = (52, -^2, ^^2)- This 



equivalent system can be described in terms of the series 
product Gt = G2 < Gi [M] , defined by 

G2 < Gi = {8281, L2 + 82LUH1 +H2+ Iir{lIS2Li]}). 

(15) 
Note in equation ( 15 ) the order of the operators is impor- 



tant. The series product provides the three parameters 
for the combined or total open system G in terms of the 
parameters for each of the systems Gi and G2. 



Go 



Gi 



G = G2 < Gi 



FIG. 4. Two quantum systems cascaded, so that the output 
of system Gi becomes the input of system G2. Importantly 
the flow of information is directional. In (b) the circuit has 
been simplified using quantum network theory to an equiva- 
lent system G = G2 < Gi. This type of network topology is 
called a cascade or series connection. 



III. SINGLE PHOTON FIELDS 

The master equation for a Markovian coupling of a 
system to a boson field in a continuous-mode one or two 
photon state was first treated in [?T]. In this section, 
we review the problem of determining the associated fil- 
ter (stochastic master equation) for an arbitrary system 
G = {8,L,H) driven by a single photon field [12]. In 
section fill G| wc generalize the master equation and filter 
to include any combination of single photon and vacuum 
as a probe field. The final section, section |IIIH[ is an 
explicit example of the homodyne single photon filtering 
equations for a two level atom. 



A. Continuous-Mode Single Photon States 

There are many ways to generate single photon states 
j47j . One common technique for creating heralded single 
photon states is by spontaneous parametric downconver- 
sion (SPDC). The photons from such a process are inher- 
ently multimodal [43], and spectral filtering is typically 
performed to get a single mode photon. 

The creation operator for a photon with one-particle 
state ^ is 



BHO 



Ht)dBHt) 



(16) 



normalized so that ||^|1 = J^ \S, {t)\^ 
photon state is then defined to be 



dt = 1. The single 



(17) 



One may interpret this is the frequency domain as |1^) = 
J^^£.{ojW{uj)\0) where ^ is the Fourier transform of £_ 

and b (w) the formal transform of the input process. This 
representation is often referred to as the multimode, or 
continuous-mode, single photon state, see for instance 
[B:, Sec. 6.3], [13 Sec. 14.2], [46, Eq. (9)]. 

Much of the calculations that follow will involve the 
identities 



dBit)m) = atmdt, 

dAit)\U)=mdBHt)\0), 



(18) 



and this will be the origin of the departure of the master 
and filter equations from the vacuum case. 



B. Single Photon Master Equation 

Without loss of generality we fix the initial state of the 
system to be a pure state po — \v){v\ ^^'^ our aim is to 
obtain a differential equation for the expectation 

for arbitrary system operator X. Starting from the 
Heisenberg-Langevin equation as before, but now using 
the identities ( 18 ) we find 



dt 



njl'{X)=En[3t{CGX))] 

+ Eoi[jt{8^X,L])]e(t)+E,o[jt{[L\X]8)]at) 

+ Eoo\jt{8^X8-X)]\^{t)\^ 

= r^l\CGX)+w'^H8^[X,L])C{t) 

{p,X]8)m+^^'iS^XS-X)\at)\' 



■ VJ. 



10 



where 



with 



E,k[A]^{Tjc^j\A\r^cl,k} 
u^l'{X)^E,k[jt{X)] 



|0), 3 

lie), J 



Rather than finding a single master equation as in the 
vacuum case, we end up with a system of equations 

ml\X) = ml\CX)+m",\8^[X,L])e{t) 

+wl%[L\X]S)m+^TiS^X8-X)\m\^ 
vul^iX) = zul'{£X)+zu^'{8^[X,L])e{t), 
va°\X) = vu^'{CX)+vu^"{[L\X]8)at), 
w^°{X)^m',\CX), 

with initial conditions 



(19) 



w^'iX) = <'(X) = (,/,Xry), ml^iX) = u^'^^X) = 



(20) 



The main feature here is that the differential equation for 
expectations tu^'^ depends on lower order w^^ , allowing 
us to solve for m^^ inductively. Likewise, defining the 
traceclass operators g''' via 

ir{g^\t)^X]^vjl''{X), (21) 

we obtain a system of equations 

ei"(t) = £*eio(t) + [V°w,i^]ew, 

eOi(t) = £*e"i(t) + [L,g"0(f)5t]r(t), 
lP\t)=Cd>\t), (22) 

with 

gii(O) = gOO(O) = |ry)(r?|, ^^"(O) = e°i(0) = 0. 

Note Q^'^it)'^ = Q''3{t). 



which means that at some stage the atom decays into 
its ground state | J,) creating a single photon in the out- 
put. The mechanism for producing the single photon is 
therefore spontaneous emission due to the coupling to 
the vacuum fluctuations. Here U- is the lowering oper- 
ator from the upper state | t) to the ground state | |). 
The Schrodinger equation for IV"*) = V"(i)| t> ® |0) then 



becomes d\tpt) 



-dt 



iV't), 



X{t)a-dBl - \\\{t)\a+u. 

and it is an elementary calculation to see that this has 
the exact solution 



iv-t) = ^/W)\ t> 



i)®Bl{0\Q) (24) 



where B^ (^) = J„ S,sdB*, and (to preserve normaliza- 
tion) w (t) = J^ 1^ (s) I ds with the complex- valued func- 
tion ^ (•) related to A (•) by 



A(t) 



V^Ht) 



e(i) 



(25) 



An Input-Output Model of Single Photon 
Signal Generation 



In section [III E| wc will set up a general technique for 
deriving the filtering equations for situations including 
the single photon input field. It is possible to give an 
alternate derivation in this case motivated by the idea 
of using a pre-interaction preparation where a vacuum 
input is first passed through a fixed system in order to 
generate the one photon field. Our motivation for consid- 
ering such a scenario stems from statistical and engineer- 
ing modelling where it is common practice to use 'signal 
generating filters' j38| driven by white noise to represent 
colored noise. Analogously, in this section, we construct 
a quantum signal generating filter M = {Sm, Lm, Hm)- 
Cascading the single photon generating filter M with the 
quantum system G we wish to probe, Figure [Sj we cre- 
ate an extended system. Because this extended system 
Gt = G < M is driven by vacuum, the master equation 
and quantum filter follow from the known vacuum case 
upon substitution of the parameters for the cascade sys- 
tem (Section [mEl). 



white 

noise 



S(() 



B„,„(f) y{t) 
>4hd1 



measurement 
detector signal 



FIG. 5. An ancilla system M is used to model the effect of 
the single photon state for B{t) on the system G. 

The idea behind the signal generating filter M is sim- 
ple. We take the filter to be a two level atom initially 
prepared in its excited state | f) . The interaction with 
the vacuum input is taken to be 



Since w (0) = ||^|| = 1, we therefore generate the limit 
state 



l^oo) = ii)®i3^(0lo) = i;)®|i 



a- 



If we wish to engineer a desired state ^, then we need 
to choose the (time-dependent) coupling strength A ac- 
cording to ( 25 ) . 



D. The Extended System 

We now define our extended system as the cascade 
system Gt = G <\ M, as in Figure [5] where using the 
cascade connection formalism from Section [II Gl we have 







{SM,LM,HM)^{I,X(t)a-,0), 



(23) 



(26) 



Let us denote by U{t) the unitary for the extended 
system driven by vacuum for the parameters Gt on the 
ancilla-|-system Hilbert space. Specifying an initial state 
I t) ^ \r]) (E) |0), we consider the expectation 

wtiA (g) A) = E^^o[U\t){A (E) X)U{t)], (27) 

(here A is an ancilla operator, and A is a system opera- 
tor). 

In order to be useful, the extended system Gt (driven 
by vacuum) must be capable of capturing expectations of 
X{t), for arbitrary operator A of the system G, at time 
t as if it were driven by the single photon field. That is, 
we must have 

E^^[X{t)]^E^„o[U\t)(I<E>X<E>r)U{t)], (28) 



that is we have the situation outUned in equation (II]) with 



with 



We are required to show that 

^^^[X{t)]=¥.^no[U\t){I®X)U{t)] 
holds for any operator X of the system G. 



(29) 



Our verification of ( 29 ) is to compare the differentials 



of both sides. Now the left hand side of ([29| is just 
the single photon expectation vjI^{X) = Eii|X(i)] = 
E^^[X(t)], whose differential equation is determined from 
the system (19). The differential of the right hand side 



of ( 29 ) may be found using the Lindblad superoperator 
Cg^\A ® X] for the extended system, which may be ex- 
pressed in the form 

£gt [A®X]=A® CgX + {Vl^,A) ® X 

+ l\jA ® S^ [X, L] + ALm ® [L"^ , X]S 
+ LlfALM<E){S^XS-X), 

for any ancilla operator A and system operator X. We 
first observe that 



Vl,Ai)^o, T^LA<y-) 



ml 

2w(t) 



(T_, 



where w(i) = /^ |^(s)p(is. Then we have 

jE^,o[UHt){I^X)U{t)]^ 

wl\CGX) + m",\S^X,L])C{t) 
+^l°{[L\X]S)m + ^^"iS^XS-Xmt)\\ (30) 



where 



~ jk 



^rix) = 



Wjk{t) 



(31) 






Qqo Qoi 

QiQ Qii 

■wqq wai 

Win wii 









Notice that equation (30) for m}^(X) has the same form 
as the m}^{X) equation in (19). In general, the equa- 
tions for m-j. {X) have the same form as equations ( 19 ) for 



nj-j. {X). Since at time i = we have TOq (X) = zuq {X), 

it follows that mj (X) = ^l [X] for all t. This estab- 
lishes the identity ( 29 ) . 



E. Single Photon Stochastic Master Equation 
(Filter) for Quadrature Phase Measurements 

In this section we explain how the quantum filter for 
the conditional expectation 

^,"(X)=E,5[X(t)|y(s),0<s<t] 

for the system G driven by a single photon field may now 
be obtained from the quantum filter for the conditional 
expectation 

^t{A®X)^¥.^,^Q[tjHt){A®X)U{t)\I®Y{s),0<s<t] 

(32) 
for the extended system Gt = G <\M driven by vacuum. 
Indeed, we have 

dMA®X) = nt{CGT{A(E)X))dt 

+{nt{A (X) XLt + L%A ® X) 
-MLt + Ll)nt(A (E) X))dW{t), (33) 

where dW{t) = dY{t) - nt[LT + L\.)dt. If we define 






Wjk{t) 



(34) 



where Qjk and Wjk{t) were defined in the previous sec- 
tion, we obtain the coupled system of nonlinear stochastic 
differential equations 



d7rl\X) ={7rl'iCX) + 7r^\SHX,L])e{t) + 7rl%[L\X]S)m + n^^'iS^XS - XM^dt 

+ {nl'iXL + L^X) + ^'l\S^X)e{t) + Tr]\XS)m - Trl\X)Kt]dW{t), 
d4\X) ^{n]\CX)+^r{S^[X,L])C{t)]dt + {{^}\XL + L^X)+n^,\S^X)C{t)^^l\X)Kt]dW(t), 
dTT°\X) ^{7T°\CX) +Tr°''{S^X,L])C{t)}dt+ {{Tr°\XL + L^X) +TT^'>iS^X)C{t) -Tr^t\X)Kt}dWit), 
d^°°(X) =Tr°'>{CX)dt + {7T°"{XL + L^X) - Tr^"{X)Kt}dWit). (35) 



Here, 



and the innovations process W(t) (given above) may be 



Kt = nl\L + it) + n^\S)m + vr,i"(5^)r (i) (36) 



expressed as 



For the extended system we have 



dW{t) = dY{t) ~ Ktdt. 



(37) 



We have 7rOi(X) 



^10 



are TTn 



{X^)\ and the initial conditions 
{X) = ^l\X) = {^,Xrj), 7riO(X) = 7r«i(X) = 0. 
In order to see that the single photon quantum fil- 
ter is given by the system of coupled equations (351, we 



must show that the conditional expectation for the sys- 
tem driven by the single photon field is given by 

TTtiX) =Er,^[X {t)\Y{s), 0<S<t] 

= E^r,a[UHt){A (g) X)U{t)\I (g) r(s), <s<t] 
^nl\X). (38) 

To obtain the filter, we again apply the characteristic 
function technique, setting C = Cg{t) as before with 
dcg{t) = g{t)cg{t)dY{t). We need to verify that 



E^o[jt{X)cg{t)] = E^o[MX)cg{t)] 



(39) 



E^r,o[UHt){A®X)Uit)cgit)]=E,j„[ntiA®X)cgit)], 

(40) 
for all functions g, arbitrary ancilla, and system operators 
A and X respectively. Hence ( 39 1 will follow provided we 
can show that 



E,k[X{t)Cg{t)] 



E,^o[U^t){Q,k®X)U{t)cg{t)] 

Wjkit) 



(41) 



However, equation (41) may be verified in exactly the 
same way we proved that wj (X) = w^ (X) in the pre- 
vious section, that is, by comparing the differentials of 
both sides of ( [4l] ). The details of this calculation are 
omitted. 

Now, write 7r^''(X) = tr{{gi''{t)YX). Then from the 
differential equations for n^ {X) and the definition g-''^{t) 
we immediate get the differential equations for the evo- 
lution of g-'''{t), as follows: 



dp^'it) = {/:>"(t) + [Sp'\t),L^]m + [L,pi°(t)5t]r(t) + (5p"°(f)^t _ p00(^))|^(^)|2|^^ 

+ {Lp'\t) + p'\t)L^ + pi°(t)5tr (i) + Sp"'{t)at) - Ktp'\t)}dW{t), 
dp'"{t) = {{C*p'"{t) + [Sp'"{t),L^]^{t)}dt + {Lp'"{t) + p''{t)L^ + Sp^'^mit) - Ktp''{t)}dW{t), 
dp'Ht) = {/:V°^(t) + [L,p'^{t)S^nt)}dt + {Lp^\t) + p°\t)L^ + p™(t)5tr(i) - Ktp^\t)]dW{t), 
dp'^'^it) = C*p'^°{t)dt + {Lp°%t) + p'^°{t)L^ - Xtp™(i)}dT4^(i), 



(42) 



where 

Kt = tT{{L + L^)p'\t)} 

+ tr{Sp"\m{t)+tr{S^p'°{t)}C{t), 

with the initial condition 

p"(0)=p°O(0) = |r/)(r;|, piO(0)=p"i(0) = 0. 



F. Single Photon Stochastic Master Equation 
(Filter) for Photon Counting Measurements 



In this section we briefiy derive the filtering equations 
for photon counting measurements. The quantum filter 
for the photon counting case is given by the system of 
equations 



dnl'iX) = {nl\CX) + 7T^\S^[X, L])Cit) + n^^ {[L\ X]S)m + n'i^iS^XS ~ XMmdt 

-f jz^r^ (nlHL^XL) + n'^\S^XL)C{t) + nl'{L^XS)m + tt^" {S^ X S)\m\') " ^t"(^)}dA^(t), 

dirl^iX) = {nl'{CX)+n^\SHX,L])C{t)}dt+!^,.^' (^^"(Ltxi) + ^™(5txL)r(^)) - ^f (X)|diV(t), 
d7:^\X) = {4\CX)+7r^'{[L^,X]S)m}dt+i^,^r' ^^(itXL) + ^,"°(Ltx5)e(t)) - ^f (X)|diV(t), 
d7r°°{X) = n°,'{CX)dt + \,.i'{n°,'{L^XL)) - 4\X)]dN{t), 



where 



u,^nl\L^L)+7:°\S^L)Cit) 
+ <(Lt5)e(i)+<"(/)|e(t)P, 



and the innovations process N(t) is given by 
dN{t) ^ dY{t) - vtdt. 
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G. Combination of One Photon and Vacuum States 

In this section we take the state of the field to be in a 
state defined by the density operator 



Pficld 



jk 



7fejl0j)(0fe| 



(43) 



where we use the notation introduced above for the pho- 
ton \(pi) = |1^) and vacuum |0o) = |0) states. The coef- 
ficients 7jfc must of course satisfy the condition that the 
2x2 complex matrix 






7ii 710 
700 



(44) 



is a density matrix, i.e. pa > 0, Tr[/9a] = 1. By choos- 
ing the coefficients jjk appropriately we can model an 
input field that is any combination of single photon and 
vacuum. For example: the single photon field is given 
by 7ii — 1 and all other coefficients are zero, a superpo- 



sition like 



"i|lf 



/ r ~ Lti |x^/ -I- aQ |0) is obtained by setting 

7ii = l"iP: 7io = ai^o' 7oi = "oai, 7oo = koP; and 
a simple combination is pficid = V |1)(1| + {^ — v) |0)(0| 
where 711 = ry, 700 = 1 - 7/ and 710 = 701 = 0. 



1. The Master Equation 

The expectation vjt{X) = {X{t)) of the system opera- 
tor X{t) when the system and field are initialized in the 
state |77)(??| fX" PfioU is given by 

Y,^,tM,k[X{t)] 
Jk 



Wt{X)=Er, 



[xm 



Il7jfei 
]k 



4'ix), 



(45) 



where wj {X) are defined in section HIIBI While there 
is no differential equation for wt{X), it can be computed 
from the weighted sum ( 45 1 , Figure pj From equation 
( 45 1 we see that the density operator for the expectation 



(46) 



TOt(X) = tv{g{t)X} is given by 

git)^J2^kjg'''{t), 

jk 

where the g'^it) are the density operators introduced in 
section IIIIBI 



2. The Stochastic Master Equation 

Turning now to the problem of determining the filter, 
we again make use of the cascade extended system from 
Section lill Dl Now we have 

■^tiQjk 'S)X) 



]k 



ix) 



Ojkit) 



(47) 







F 


F 


r 


r 



> Wt{X) 



master 
equations 



weighting 



FIG. 6. The expectation ■cut{X) = {X{t)) of the system op- 
erator X{t) = jt{X) when the system and field are initialized 
in the state \rj){ri\ (g) pfieid may be calculated by weighting the 
solutions zul (X) from the single photon master equations 
(191). 



and so if we define the matrix 



i?w = E 



jk 



Wjk{t) 



ijk, 



(48) 



where Wjk{t) and Qjk are as defined in Section 



have (using (l47|), (|45|) and (l27|)) 

wt{X) ^ wt{R{t) ® X) 



HID 



we 



(49) 



Note that the definition (27) oiwt{R{t)®X) involves the 
ancilla system initialized in the excited state |e) = | f)- 
The conditional expectation 



7r,(X)-E,,.,^„,[X(t)|r(.),0<.s<i] 



(50) 



corresponding to the field in the state Pfidd is related to 
the conditional expectation TTt{A (S) X) for the extended 
system (see (32)) by the Bayes relation 



MX) 



nt{R{t)(E)X) 



(51) 



ntiRit)(g>I) 

Division by the denominator in (|51|) is needed to ensure 
the normalization 7r((/) 
show that TTt{R{t) ® X) 
lently 



1. To prove (51), we need to 
TTt{R{t) ® /)7rtX^), or equiva- 



Et„o[^t(i?(t) ^ X)cg{t)] = E^,,oMR{t) ® I)MX)cg(t)] 

for all choice of characteristic functions Cg (i). How- 
ever, E^jjo[MR{t) (S) X)cg{t)] equals E^^o[UHt){R{t) ® 
X (g) T)U{t)cg{t)], but by the extended system represen- 
tation this is just E^pj.^jj[X(t)cg(i)] which in turn equals 
E^p^[ntiX)cgit)] = Ello[UHt){R{t) ® I)U{t)TTt{X)cg{t)] 
which establishes the Bayes relation (51). 
Since 



MRit)®X)=Y,^,,7rt{X), 

jk 



(52) 
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Jk 



where irj {X) is defined by (34), the desired conditional 



expectation may be expressed as 



MX) 



Y.i.kljkTTl (!) 



(53) 



Again, there is no differential equation for ■nt{X)\ in- 
stead it is computed from a normalized weighted sum, 
(53), and the filtering equations (35), Figure [t] The cor- 
responding conditional density operator is given by 



P(0 



E,fc7fc,tr{p^'=(f)}' 



(54) 



where the conditional quantities />''^(i) may be computed 
from the single photon filtering equations ( 42 ) . 



Y{t) 



measurement 
signal 





^"(^) 








w 








w 







-> 7rt(X) 



filtering 
equations 



weighting and 
normalization 



FIG. 7. Relationship between the measured signal and the 
filtered estimate. The differential equations ( 35 1 must be in- 



tegrated to compute vr^ {X). Then depending on the state of 
the input field, the probe, these estimates must be weighted 
by t he appropriate coefficients and normalized as specified in 
(53 1 to produce the desired conditional expectation 7rt(X). 



These expressions allow filtering on any combination 
of a single photon and a vacuum state. One notable case 
it that of simple combination of one photon and vacuum 
(Pprobc =p|l)(l| + (1 — p) |0)(0|) which is an experimen- 
tally accurate model for the output of the SPDC process 

HI- 



1. The Master Equation 



As there are only two components to the input field 
1^) we only need to consider a generalized Bloch repre- 
sentation of the form 



rjk 



(^(Jkj _^ ^Jk^^ _^ yJk^^ _^ ^jk^^^ 



(55) 



with j,k e {0,1}. Note that x°", y"°, z°" and x", y 



nil ,,11 



vii 



are real, while a; , y % z may be complex. For the 
master equation the coefficients cP'' = 1, we will however 
require them later. Also note, for example, tr {g'^'^Cx} = 



^00 



^01* 



and tT{g"^'^ax} 
Using equations ([22|) derived in section IIIB 



compute the generalized Bloch component equations for 
the unconditioned evolution of our system. We obtain 
nine coupled equations for the nine components: 



^00 



-2ujy 



00 _ „oo 



K 



yO" = 2^x"°-^y°", 



i00^„^(1^^00) 

^01 _ '^^oi o, ,„,oi ^Tt('+^*~oo 



a;"^ ^ __a;Oi _ 2^yOi _ ^^^tyz°'', 

f'=2u;x''~^y''~^V^m*^'', 
iOi ^ _^^oi _ , 



'm* + ^v^y'"atr 



i" = -^a;" - 2c^y" + V^z°'^{t) + VKz°'*^{tr, 

y" = 2ujx^^ - 1^" + tv^z'^^^it) - t^z°^*atr, 

i" = -K - Atz" - y/KX°^^{t) - iy/Ky°^^{t) 



These equations must combined using ( 55 ) to calculate 
relevant quantities. 



H. Illustrative Example of Single Photon Master 
and Filtering Equations 

Here we give an explicit example of the filtering 
method derived above to a system driven by a pure single 
photon with homodyne measurement of the amplitude 
quadrature of the field. We take the system of interest 
to be a two-level system (a qubit). In particular, we con- 
sider a damped qubit L — ^fka- with internal dynamics 
given by i? = wCTz and no scattering i.e. S — I. Here 
K > is a scalar parameter often referred to as the mea- 
surement strength and a; is a frequency. 



2. The Stochastic Master Equation 



Now we turn our attention to calculating the filter- 
ing equations for this example. Recall that c^'' = 1 in 
the generalized Bloch representation for the master equa- 
tion. For the filter, the coefficients c^'' are not necessary 
unity, though we always have c^^ = 1 at all times. How- 
ever, unlike the master equation, this will not be so for 
c°^, c^°, c°°, as these coefficients will now evolve in time. 
Now, substituting our particular choice of {S,L,H) into 
equations ( 35 ) , the quantum filter for the two-level sys- 



tem is given by the finite set of coupled equations 



12 



dx 



" ={ - fa;" - 2c.y" + ^^z^'m + ^nz^'*m*] 



di 



rfy" ={ 



2w2:^ 



-icOiew + icOi*ew*)a;" 



-y" + i\/)^z"i^(t) -i^)^. 



KZ 



'i*^(t)*} 



dt 



y'^'*m*+y''m-{v^^' 



1^01 



eW 



1„01* 
2^ 



ew*)y 
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dVF(i), 



^11 



dz" = 



Kxioe(i) - i\/^y^°e(0 






1„01 



ew 



icOi*CW*)^" 



df 



dW{t). 



IV. SUPERPOSITION OF COHERENT FIELD 

STATES 

In this section we turn to the problem of determin- 
ing the master equation and the quantum filter for sys- 
tems drive by a boson field whose state is a superposition 
of continuous-mode coherent states. In section IIV Al we 
describe continuous-mode coherent states and superpo- 
sitions of them, as well as the action of the quantum 
noises on such states. Section IIVBI is devoted to the 
derivation of the master equation for superpositions of 
coherent states. In section |IVC| we develop a cascaded 
system signal model. This model allows us to use the 
methodology from Section [lll[ with appropriate changes 
due to the nature of the superposition of coherent states, 
to derive the filtering equations in section |IVD| Then 
we give the filter for the case of photon counting in sec- 



tion IV E and generalize to mixed input states. 



dc' 



,01 _ I y^x°i+c™^(i)*-(V«^" \,^(.. 



dx 



°' = { -f a;Oi - 2wyOi - V^^(i)*zO" } 



dt 



nil 



dy°^ = [2iux°^ - f yOi _ ^^^(^ty z°°} dt 
+ 1 2/°"^W* - (v^-"+ ]dW(t) 
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dc™ 
dx"o 






^x°° - (\/«2;" I ^|y/^^ 



{-2..y°"-V°} 



dt 






dW{t), 



dy"0 = {2L.x™--y™}di 

The innovations process is given by 

+c"^it)at) + c^%t)Cit))dt. (56) 



A. Superpositions and Combinations of Coherent 
States 

Typically single mode coherent states of a field are de- 
noted by I a). In this paper we shall often refer to a 
superposition of continuous-mode coherent states as a 
(continuous- mode) cat-state [HI HHj- Formally, the su- 
perpositions of continuous-mode coherent states is given 
by 



i^> = E 



Si a. 



(57) 



3 = i 



where \aj) are coherent states, determined by func- 
tions aj{t) with aj ^ ak if j ^ k. The superposi- 
tion weights Sj are complex numbers such that (?A|V') — 
Si k ^*j^k{<^j\<^k) — 1 (i-e., i/" is normalized and is a pure 
state vector of the field). Given a function a, the co- 
herent state I a) of a continuous- mode field is given by 
the displacement or Weyl operator D{a.) applied to the 
vacuum state of the continuous-mode field: 



D{am. 



(58) 



The inner product of two coherent states |a), |/3) in the 
Fock space is given by 



a, 



where 



exp 



1 



1 



||/3|r+(a|/3)), (59) 



is the LP' norm. The normalization condi- 



tion for the superposition state (57) means that the co- 
efficients must satisfy '^j k^*j^k9jk = 1, where gjk = 
{aj\ak). 

More generally, we may consider a field density opera- 



tor 



Pfioid ^^lkj\aj){ak\ 

jk 



(60) 
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that generalizes the superposition state \ip) to allow for 
statistical combinations of coherent states. The normal- 
ization for the state pfieid is J2j.k "fjk9]k = 1- 

In what follows the action of the quantum noises dB 
and dA on coherent states will be important: 



dB{t)\a) =a{t)\a)dt, 
dA{t)\a{t)) =dB*{t)a{t)\a). 



(61) 



and g^'^iO) = \v){v\- 



The master equations (65 1 and (68) consist of a 



weighted sum of cross-expectations. Clearly these equa- 
tions reduce to the vacuum master equation if the only 
term in the superposition or combination is the vacuum. 



C. Extended System 



B. Master Equation for Systems Driven by a Field 

in a Combination or Superposition of Coherent 

States 

Again, before we derive the master equation we intro- 
duce some notation that helps to formulate the master 
equation. Recall that we defined the asymmetric expec- 
tation Ejk [X IS) F] = {r]\X\r]){(j)j\F\4>k). In section["^ 



III 



took the field states \(j)j) , \(j}k) to be either vacuum or one 
photon. In this section we use this same notation but the 
field states are understood to be continuous-mode coher- 
ent states i.e. \aj),\ak). The indices j,k now take the 
values 1, . . . ,n. 

The expectation of an arbitrary system observable, 
with respect to the state |?7)(77| (g) padd, at time t is 

mt{X)^E,p,^jX{t)]. (62) 

Using the notation (similar to the single photon case) 

wt{X) = Ejfc [Xit)] = {r^a,\Xit)\r^ak), (63) 



with Pficid as given in |60[ we may write ( 62 ) as 



■^t{X) =^7jfctA7^''(^)- 
jk 



(64) 



As in section [III B[ we can derive the Heisenberg mas- 
ter equation by taking the expectation of the equation of 
motion for an arbitrary system operator dX{t), i.e. (|4|. 
Doing so yields the equations 



ujr(x)=ujr{grx), 

where we define a new superoperator 

g^X = CX + S^[X,L]a*it) + [L^,X]Sak{t) 
+ {S^XS-X)a*{t)ak{t), 



(65) 



(66) 



with initial conditions zuq (X) — {ri\X\r])gjk. Note that 
equations ([65|) are uncoupled. 



The corresponding density operator is 

B(t)^Y.^,kg'Ht) 
Jk 

where 



(67) 



g^'^ = gt'[g] = C*g+[Sg,L^]a,{t) + [L,gS]al{t) 

+{SgS^-g)a,it)alit), (68) 



In this section we describe a cascade extended system 
Gt — G <\ M that will be used in section |IVD| to deter- 
mine the quantum filtering equations for the mixed or 
superposition of coherent state field. The ancilla system 
M will be an n-level system, with orthonormal basis |j), 
j =^ 1, . . . ,n. The parameters for this system are 



where 



M^{I,Lm,0), 



^M = 5I«jWli)01. 



(69) 



and we take the initial state of the ancilla to be the den- 
sity matrix 



""" Jk 



Pa = 



(70) 



where Na — X); In ^^ ^ normalization factor. The ex- 
tended system is 

Gt = G<M :^{S,L + SLm,H + Im {L^^SLm})- 

Define Qjk — \j){k\- Then a straightforward calcula- 
tion shows that 



CL^iQjk) =mjk{t)Qjk, 



(71) 



where 



m,k{t) = a*{t)ak{t) - \\a,{t)\^ - ^\ak{t)\\ (72) 

Now consider the extended system Gt initialized in the 
state Pa <E)\r]){ri\(E)\0){0\ (driven by vacuum |0)). Then the 



methods used in Sections |III D| and |III G| may be adapted 
;ase to si 



to the present case to show that 

■^tiQjk <E)X 



Wjk{t) 

where Wjk(t) is defined to be the solution of 

Wjkit) = mjk{t)wjk{t), Wjk{0) = 



Na9jk ' 



(73) 



(74) 



and 



E^,,^jXit)] ^E,^,o[uHt)iRit) ® X)Uit)], (75) 
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where 



R{t) 






(76) 



These expressions are very similar to the photon case, 
but with some important differences. For instance, the 
ancilla was initiahzed in the excited state for the photon 
case, while here for the mixed coherent case the initial 
ancilla state is the density pa- 



D. The Stochastic Master Equation (Filter) for 
Amplitude Quadrature Measurements 

The quantum filter for the general combination of co- 
herent state case may now be derived in exactly the same 
way as was done for the combination of single photon and 
vacuum in Section [Til G[ The conditional expectation we 
are interested in is 



TTtiX) =E,,p,^jX{t)\Y{s),0 < s <t], 



(77) 



where now pgcW is given by (60). Equations (51), (53), 



and the new superoperator Jij (•) is defined by 

Hf (X) 



nl''{XiL + Sakit)) + iL^+S^a*it))X) 
111 

Na 



w'(^) E ^'^"(^ +^^ + ^"' w + s^^ii*))- 



As before, we may write tt^ (X) = tr {g^''{tyx^, where 
0'^{t) satisfies the coupled differential equations (for 
j,k = l,2,...,n): 

dp'^'it) = gf '^[/^■'=(i)]di + Hf * [p'^'it)] dW{t) (78) 

where 

nt [fyi^] ee(L + Sau{t))p^'' + p^^L^ + S^o^*{t)) 



v" Y.w^' [(^ +^^ + ^"'(^) + ^^"' (^))/'"] 



with initial conditions p^ (t) — \Ti){ri\gjk (recall that 
9jk = ( aj \ak))- The conditional density o pera tor is given 



and (154|) again hold, but with modifications to the terms by (IM), with the p'''{t) given instead by ((tsI 



as described above. The filtering equations are as follows. 
The conditional quantities ttj (X) satisfy the coupled 
system of equations 

dntix) = nt{g^^x)dt + nl^ {X) dw{t) 

where the innovations process W{t) is a Wiener process 
and is given by 

dW{t) ^dY{t) 



We remark that the innovations for the cat case now 
depends on the weights, in contrast to the mixed pho- 
ton/vacuum case. 

E. The Stochastic Master Equation (Filter) for 
Photon Counting Measurements 

Analogously, we may also compute the quantum filter- 
ing equations for a system driven by a coherent super- 
position in the case where the measurement performed 
on the output field, Y{t), is photon counting. The filter- 
ing equations in the Heisenberg form are given by (for 
j,k = l,2,...,n): 



., .. ., /7:i''{L^XL + ak(t)L^XS + a*(t)S^XL + a*{t)ak(t)S'fXS) -. \ 

dnl'^iX) = ^f (g-'^(A))dt + -^^ ^ '— jAU^ L _ ^Jfe(x) ]dN(t), 



where 



dN{t) = dY{t) - Y^ ^Trf (i^L + aj{t)L^S + a*{t)S^L + \aj{t)\^I)dt, 



j=i 



(79) 



jkf 



and with initial conditions ttq (X) — {ri\X\ri)gjk 



F. Illustrative Example of Superposition of 
Coherent State Master and Filtering Equations 



r 



wave field to be in a Schrodinger cat state of the form 



I*) = _^(ci|ai)+C2|a2)) 



(80) 



We give an explicit example of the filtering method 

derived above to a system driven by a continuous-mode where ci, C2 are complex numbers that determine the 

cat state under homodyne measurement of the amplitude relative weighting and phase between the two complex 

quadrature of the field. We consider the input travelling (continuous- mode) coherent state amplitudes ai and 0^2, 
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and Af is real normalization factor. We take the system of 
interest to be a two- level system (a qubit). In particular, 
we consider a damped qubit L — ^/kcf- with internal 
dynamics given hy H = loUz and no scattering i.e. S = I. 
As there are only two components to the input field 
I ^) , we only need to consider a generalized Bloch repre- 
sentation of the form g^'^ = {c>^I+x^^crx+y-'^(Tji + z^^az) 
with j, fc e {1,2}. Using equations (68 1 and (67 1 derived 



earlier in this section, we now compute the equations for 
the unconditioned evolution of our system i.e. the master 
equation(s). Thus the master equation is given by the set 
of equations: 



^11 



•,11 



i" 



Kx" -2wy" + y;^z"[ai(t)+a*(t)], 

-^Ky" + 2wx" -I- i^fflz^\ai{t) - al{t)], 
-k(1 + z") - ^^a:"[al(i) + al{t)] 



Then the propagating equations are: 
dc" =(\//^a;" + {ai{t) + al{t) - i/)c")dW^, 
dx" = (-^Kx"-2a;2/"+v^z"[ai(t)-t-ai(i)] jdi 
+ (\//^(c" + z") -I- (ai(t) + a* (t) - i^)x^^)dW, 

dy^^ = (--Ky" + 2tja;"+iv^z"[ai(i)-ai(t)] j dt 

+ {ai{t) + al{t) - iy)y^'^dW, 
dz^^ ={-k(c" + z") - ^/^x"[al(t) + a\{t)] 
- i0«y"[ai(t) - al{t)\}dt 
+ (-x" -f- (ai(i) + a*(t) - v)z^^)dW- 

dc^^ ^{x^'^ + (a2(i) + al{t) - i^y^)dW, 
dx^^ = [-\kx^^ - 2a;|/^2 + ^/Kz^^[ai{t) + a^(i)] j dt 
+ {y/^{c^^ + z^^) + (a2(t) + a\{t) - i^)x^^)dW, 



±12 = 



y'' = 



±21 



f' = 



--Ky^2 + 2ujx^^ + iy/^z^^[ai{t) ~ a;(t)], 

--Kx^i - 2wy2i + y^z2i[a2(t) +a*(t)], 

ny^^ + 2ujx^'^ 

-^-i^/Tiz [a2{t) ~ a\{t)], 

+ai(t)] - i^/Ky'^^[a2{t) - ai(i)], 



±22 = --^22 - 2cjy22 + ^z^'^[a2{t) + al{t% 
y22 = -^ny^^ + 2cja;22 _^ i^222[a2(i) - ^^(t)]. 



i22 



k(1 + z22) - y^a;22[a2(t) -f a2(t)] 



-~i^y^^[a2{t) - al{t)l 

and the initial conditions are s^^{Q) = trjcTs |?7)(?7|} for 
s G {x, y, z}. It is easy to see that the coefficients c'^{t) = 
1 at all times by making the substitution X = / in ( |66[ ) . 
Next we present the coupled conditional equations which 
propagate the filter. Let 



dyi2 ^ — ^yi2 ^ 2ujx^'^ + i^/llz^'^[ai{t) - a*2{t)] ] dt 
+ {a2{t) + al{t)-v)y^^dW, 

dz^^ = {-k(c12 + zi2) _ ^x^^[ai{t) -t- alit)] 
- i^/Ky^^[a2(t) - al{t)]}dt 
+ (-a;i2 + (Q,2(t) -I- a*(i) - i.)zi2)diy; 



dc^i =(v^x2i -I- (ai(t) -K a*2{t) - vY^)dW, 

= (--Kx^i ~2w2/2i +y^z2i[ai(t)-|-a;(i)]) dt 
+ (\/^(c2i -I- z2i) -I- (ai(t) + a*2(t) - v)x'^)dW, 
( --Ky2i _^2wa;^i + «VKz^^[a2(t) - ai(t)] j dt 

/ -k(c2i -f z21) - v^a;2i[a2(i) -f ai(i)] ^ 
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dx 



dy''^ 



dz'^^ 



-iy^y^^[a2{t)-al{t)] 
+ (-x^i + {ai(t) + a*2{t) - v)z^^)dW; 

dc"^ ^{^x^^ + (a2(t) + a;(t) - iy)c^^)dW 



dt 



|cir + |c2r 

+ |ci|2(ai(t) + at(t))c" + \c2\'{a2{t) + am^')- 



dx^'^ = I --Ka;^^ - 2wy22 + y/Kz'^'^[a2{t) + ^^(t)] | dt 
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+ {a2{t) + a;{t) - ,y)y^^dW, 

- iy/Ky^^[a2{t) - al{t)]}dt 

The innovations process is 
dW{t) = dY{t) -Y^^,[\j){j\ L + a,{t) \]){j\ ® I 

3 

+ \j){]\® L^ + a*{t)\j){j\® I)dt, (81) 

= dY{t)- [7ri(L + L^) +^aj{t) + a*(t)\dt 



2 „*, 



C,Cfc 



+ ^(a,(t)+a*(t))jdt, 



(82) 



this general combinations of the vacuum an a one pho- 
ton state through a straightforward weighting procedure. 
The filter equations themselves have potential applica- 
tions to areas such as shaping wave packet for maximal 
/ minimal absorption by, for instance, a two level atom, 
or to controlling the system so as to shape the outgoing 
field. 



We have also derived the quantum filter for cat states. 
While the concept of an environment being in a superpo- 
sition of states may seem unphysical from the perspective 
of macroscopic superselection rules, as we have seen this 
may effectively be what happens internally once a stan- 
dard input is first fed through an appropriate filter sys- 
tem M. This leads naturally to questions of decoherence 
[48] , and whether preparing input in a cat state is ad- 
vantageous in preventing decoherence of cat states for a 
given system. It is now experimentally possible to isolate 
quantum systems sufficiently well to create cat states in a 
laboratory \^-A,- The cat-state filtering equation will be 
of importance for investigating questions as to whether 
such superpositions may protected via appropriate envi- 
ronment engineering. 



and the initial conditions are s^^{Q) = tv{(Js\Tf){'q\}gjk for 
all s G {c, X, y, z}, with Uc = /. 



V. CONCLUSION 

We have shown that quantum filtering may be ex- 
tended beyond the Gaussian input situation to consider 
a range of non-classical states that are of current inter- 
est. Photon wave packet shaping is already being applied 
experimentally and our filtering equations for the single 
photon input completes the problem addressed by Gheri 
et al. in [H] by giving the quantum trajectories asso- 
ciated to the master equation they derive. We extend 
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